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ABSTRACT 

We study the scaling relations between global properties of dwarf galaxies in the 
Local Group. In addition to quantifying the correlations between pairs of variables, 
we explore the "shape" of the distribution of galaxies in log parameter space using 
standardised Principal Component Analysis (PCA), the analysis is performed first in 
the 3-D structural parameter space of stellar mass M* , internal velocity V and char- 
acteristic radius i?* (or surface brightness /.t,). It is then extended to a 4-D space that 
includes a stellar-population parameter such as metallicity Z or star formation rate 
M*. We find that the Local-Group dwarfs basically define a one-parameter "Funda- 
mental Line"(FL), primarily driven by stellar mass, Af*. A more detailed inspection 
reveals differences between the star-formation properties of dwarf irregulars (dl's) and 
dwarf ellipticals (dE's), beyond the tendency of the latter to be more massive. In 
particular, the metallicities of dl's are typically lower by a factor of 3 at a given Af, 
and they grow faster with increasing M*, showing a tighter FL in the 4-D space for 
the dE's. The structural scaling relations of dl's resemble those of the more massive 
spirals, but the dl's have lower star-formation rates for a given M^ which also grow 
faster with increasing M». On the other hand, the FL of the dE's departs from the 
fundamental plane of bigger ellipticals. While the one-parameter nature of the FL and 
the associated slopes of the scaling relations are consistent with the general predictions 
of supernova feedback (jDekel fc Wool 120031 ). the differences between the FL's of the 
dE's and the dl's remain a challenge and should serve as a guide for the secondary 
physical processes responsible for these two types. 
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1 INTRODUCTION 

1.1 Scaling Relations, Fundamental Distributions 
and Galaxy Formation in Giant Galaxies 

In the standard ACDM picture for galaxy formation, galax- 
ies form from the gas that cools within a dark matter halo 
into a disc of stars whose specific angular momentum is 
conserved (|Fall fc Efstathioulll980l ). Since the gas initially 
shares the spatial distribution and dynamical properties of 
the dark matter halo, the structural properties of the result- 
ing gaseous disc, specifically the mass of the disc of stars M, , 
the size of the disc R, and its rotation velocity Kot , are ex- 
pected to correlate strongly with those of the halo. The halo 
itself is expected to be characterised by the virial theorem, 
i.e. Mvir oc T^^ij oc -Rvir- The actual scalings of the stellar 



disc will be affected by the physics of t he baryon gas and it s 
interactions with the dark matter halo (|Dutton et al.ll2007l ). 
Thus the observed scaling relations between these structural 
parameters in galaxies help to shed light on the physical pro- 
cesses governing their formation. 

Several recent studies of scaling relations of galaxies 
have contributed to significant improvements in our ap- 
preciation of galaxy formation physics. For example, the 
scaling relation between M, and the stellar surface den- 
sity /i, (related to i?,) for 123,000 galaxies of all types in 
the Sloan Digit al Sky Survey (SPSS) sh ows a transition at 
~ 3 X 10^° M© (JKauffmann et al.|[2003al '). Above the transi- 
tion, fit, only weakly depends on M* while below it, fi, scales 
as M^'^ , consistent with si mple theoretical pr e dictions of su- 
pernova feedback (SNF) (JDekel fc Silklll986l : IPekel fc Wool 
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I2OO3I ). The transition itself can also be explained by consid- 
ering the energy co nstraints of gas ejection by supernovae 
l|Dekel fc Silklll986l ). 

The interpretation of the scaling relations between Kot 
and lumino sity L, and R* and L for ~ 1300 local spiral 
galaxies by ICourteau et alj (|2007a|) also suggests that the 
interaction of the disc and halo during galax y formation is 
a fun damental constraint to ACDM models IjDutton et alJ 
I2OO7I ). 

The study of scaling relations for elliptical galax- 
ies has re vealed the existence of a "fundamental 
plane" (FP) (iDiorgovski fc Davislll987l : iDressler et alii 19871 : 
iBernardi et al.ll2003r ). These galaxies lie on a plane in the pa- 
rameter space spanned by the logarithms of L, R, and the 
velocity dispersion a. The projections of this plane onto the 
a — L and R — L plane s are known as the Fa ber-Jackson 
(IFaber fc JacksonI 1 19761 ) a nd the Kormendv iJKormendyl 
I1977I ) relations. Recentlv. IZaritskv et all l|2006bl lal. 120071 ) 
proposed that the slight tilt of the FP can be straightened 
out essentially by adding mass-to -light ratio as an ext ra di- 
mension to the parameter space. iDekel fc Coxl (|200y) used 
hydrodynamical simulations to show how the properties of 
the fundamental plane of elliptical galaxies, including its tilt 
with respect to virial theorem predictions, can be explained 
by considering the dissipative properties of the major merg- 
ers that form them. 

Besides the structural parameters (Af, , fi,, i?*, V^), the 
physics of the gaseous systems that later become galaxies 
will also affect the global metallicity Z, star formation his- 
tory, and consequently the current star formation rate A/* 
(SFR) of these galaxies. Since Z and Af* are related to star 
formation we call these the "SF" quantities for brevity. 

Studies that involve the SF parameters with stellar 
mass have also yielded insight into the gas and stellar pop- 
ulatio n physics of galaxies. For example, iTremonti et al.l 
(I2OOJ) show from SDSS measurements of galaxies with 
lO^'^M© < M, < lO^^'^M© that Z decreases systemati- 
cally with decreasing Af, , while higher mass systems only 
show a weak dependence on Af, . The se scalings are con- 
sistent with SNF model predictions (lEiekel fc Wool l2003l: 
iTassis et al.l 120031 ). As another example. iBrinchmann et al.l 
l|2004l) sh ow trends of curr ent star formation rate Af, with 
Af», and iKauffmann et al.l (2004) show environmental de- 
pendence of this relation. This poorly explored scaling rela- 
tion may reveal insight into the effect of the disc self-gravity 
on gas processes that combine to form stars. 



1.2 Dvirarf Galaxies 

These investigations of the scaling relations and the FP have 
led to major advancements in our understanding of galaxy 
formation with respect to giant galaxies, and this motivates 
a similar investigation of dwarf galaxies. Dwarf galaxies are 
the central players in one of the most puzzling problems in 
the ACDM picture of galaxy formation, namely, the "miss- 
ing dwarf problem." This relates to the apparent discrepancy 
between the predicted distribution of halo masses at the 
low-mass end and the relat ively few dwarf galaxies actually 
observed - see for examp le IjKlvpin et al.lll999l : iMoore et al.l 
ll999l : IStoehr et al.ll2002l ). Similarly, the luminosity function 
of galaxies below the Schechter's characteristic L, is ob- 
served to be flatter than predicted in ACDM cosmology. 



implying that the ratio of ste llar mass to total mass de- 
creases with decreasing Aftot . iDekel fc Wool (2003 ) used this 
systematic variation in Af, /Aftot to predict the scaling rela- 
tions of fj,* oc Mt/R? and V with Af, for dwarf g alaxies in 
the SNF scenario. (See also lMaller fc DekelllJOOa who show 
how SNF helps solve th e angular momentum problem, and 
iDekel fc BirnboimllJOOel who detail the role of SNF in smaU 
galaxies). A dependence of the ratio of Af*/Af on M for 
low-mass galaxies, and thus a dependence of the overall star 
formation efficiency on Af , should also have consequences on 
the scaling relations of the SF quantities with Af, . 

The theoretical developments concerning the scaling 
relations between the properties of dwarf galaxies mo- 
tivate an investigation of the fundamental distribution 
of the dwarf galax ies in parameter space. For example, 
iBender et aL (Il992h plot the gal axies in what they call 



« -space. 



Prada fc BurkertI (|2002l ) (hereafter PB02) find a 



"fundamental line" (FL) of dwarf galaxies in the space of 
Aftot /I/, metallicity [Fe/H], and the central surface bright- 
ness Ev, which they explain with a simple chemical en- 
richment model. PB02 fin d that this FL is independent 
of Hubble type (see also ISimon et al.l 120061 ) . Similarly, 
IVaduvescu et al.l ((20051) find a fundamental plane of dPs in 
the parameter space of absolute K magnitude Mk, cen- 
tral surface brightness Tik and H I line width W . The 
same authors find that blue compact d warfs (BCD's) lie 
on the same fundamental plane as dP s JVaduvescu et al 



20061). Ide Carvalho fc Diorgovskil (|l992l ) and lZaritskv et al 



(2006a) find that the dE's lie in the same "fundamental 
manifold", the latter putting them on the same manifold 
as larger elliptical galaxies. However these fundamental dis- 
tributions lie in parameter spaces whose axes are various 
combinations of structural and SF quantities. In order to 
better understand the underlying physics behind these fun- 
damental distributions, a disentangling of these axes into 
more "fundamental", or linearly independent quantities is 
in order. 

The scaling relations between the structural and SF 
quantities in the dwarf galaxies were first presented as 
lumin osit y-dependent relations ; e.g.V-L: Faber fc JacksonI 



Bender et alJ 



1992 



Il976l and iTullv fc Fisheilll977l: S - L 

[Fe/H]-L, [0/HJ-L: IZaritskv et all 1 19941 and [Richer et all 
|199S . However, for the purpose of constraining and com- 
paring with physical models, the use of physical quantities 
is preferable. 

We wish to improve on these investigations. Our goals 
are two-fold: (i) present the scaling relations of dwarf galax- 
ies in physical quantities; and (ii) , find a fundamental distri- 
bution in a parameter space spanned by linearly independent 
physical quantities that can constrain physical models. We 
will show that in such a parameter space the dwarf galaxies 
lie in a linear distribution. Such a distribution, or FL, would 
be described by a one-parameter family of equations which 
narrowly constrains the possible physical scenarios of their 
origin. We will show this parameter space to include both 
structural and SF quantities. 

The dwarf galaxies of the Local Group (LG), whose 
stellar masses range from as low as 10^'^ Mq to the SDSS 
transition scale of lO'^'^'^M©, provide an excellent opportu- 
nity to achieve our goals. In this paper, we compile data 
for LG dwarf galaxies and present scaling relations for the 
structural and SF properties of the dwarf galaxies in phys- 
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ical units. We use a Principal Component Analysis (PCA) 
to search for linear distributions in structural and SF space 
and quantify their tightness. We identify the space, involving 
both structural and SF quantities, in which the LG dwarf 
galaxies form a fundamental line. 

The outline of this paper is as follows: the data base is 
described in JH followed in H2.ll by our derivation of stellar 
mass and other global properties. In 33]we describe the tools 
of our analysis, namely the fitting prescription and principal 
component analysis. In 52) we present the scaling relations 
of LG dwarf galaxies in structural parameter space. In ij5l 
we augment the scaling relations with SF parameters and 
mass to derive the fundamental linear distributions of the 
dwarf galaxies in the parameters space of structural + SF 
parameters. Our results are discussed and summarised in JHl 



2 DATA 

Our compilation of relevant global parameters for LG dwarf 
galaxies is presented in Table [l] It includes the: 

(a) stellar mass-to-light-ratio M, /L; 

(b) stellar mass M, in solar masses; 

(c) central surface density /i, in M0 kpc"'^; 

(d) exponential scale length _R* in kpc; 

(e) characteristic rotational velocity V in kms~^, taken 
to represent the halo potential well; 

(f) H I gas mass in units of lO^Mo; 

(g) mean metallicity Z; 

(h) current star formation rate M, ; 
(i) colour B — V . 

These global parameters are derived pri marily from the 
compi lations of Matoo ( 1998, hereafter M98 ) , Ivan den Berghl 
(|2000l . hereafter vdBOO), and iGrebel et al.l (|2003l . hereafter 
GGH03). GGH03 also provide absolute V magnitudes, sur- 
face brightnesses and metallicities ([Fe/H]). Of all the galax- 
ies listed in these sources, we excluded the Sagittarius dwarf 
spheroidal given its strong tidal perturbations with the 
Milky Way. 

We have classified the LG dwarfs into three basic 
types, simplifying the classification system adopted by M98, 
vdBOO, and GGH03. These are dwarf irregulars (dl), com- 
prising both dirr and Irr galaxies, but mostly the former; 
"early-type" dwarf galaxies (dE) , comprising both dwarf el- 
lipticals and dwarf spheroidals, but mostly the latter; and 
transition galaxies (Tr), normally classified as "dIrr/dSph". 
The dE galaxies tend to be less massive, less gas-rich, and 
more metal-rich than the dl's. The Tr galaxies appear to 
be a morphological transition between dirr's and dSph's. 
The Tr's tend to be currently forming stars like the dirr's, 
but are more similar to dSph's in size and shape (e.g. 
ISandage fc HoffmanI 119911 ). 

The first four parameters (M*, /x,, _R*, V) in the above 
list are the structural parameters, while Z and M* are SF 
parameters since they are related to the star formation of 
the galaxies. We describe our derivation of these parameters 
below. 



2.1 Stellar Mass Derivation 

2.1.1 Methods 

We consider two basic methods to estimate the stellar mass- 
to- light (Mt/ L) ratio of a galaxy based either on colours or 
inferred SFH's. We use a combination of both methods in 
two difi'erent ways to produce two data sets. 

The first met hod uses the cal i bration of co lours with 
model Af,/L by iBell fc de JoiiS (|200lf ) and iBeU et all 
(2003). Hereafter, we will refer to this colour- M» /I/ calibra- 
tion as "BdJ". BdJ generated model galaxies with model 
star formation histories (SFH) following a variety of star 
formation laws based on the s i mple stellar population (SSP) 
models of iBruzual fc CharlotI (|2003l , hereafter BC03). 

Our second method uses inferred SFH's rather than 
colours to calculate the M^/L ratio directly from the SSP 
models from BC03. The SSP models describe the photomet- 
ric evolution of single, instantaneous bursts of star formation 
with different metallicities. M98 provides SFH's in the form 
of histograms of relative star formation rate (SFR) as a func- 
tion of age for many of the dwarf galaxies (see his Fig. 8). 
Sampling and systematic errors in the SFH's are not easily 
estimated. However, M98 noted that for the majority of the 
galaxies, the duration and relative strength of most of the 
star formation (SF) episodes are "fairly well determined", 
while only a few galaxies (GR 8, Tucana, DDO 210, Sextans 
B, and M 32) are dominated by SF episodes whose duration 
or relative strength are "greatly uncertain" (see the caption 
of his Fig. 8). 

To calculate M^/L using SFH's from M98, metaUicities 
an d the SSP we follow the isochrone synthesis technique 
of Bruzua l fc CharlotI l| 19931 ) which assumes that a stellar 
population with a given SFH can be decomposed into a series 
of instantaneous bursts of SF. The details of our calculation 
are presented in Appendix A. 



2.1.2 Comparison of Methods 

The Mt/Lv ratios calculated with the two methods above 
are plotted against absolute V magnitude in Fig. [T] There 
is no apparent trend in either plot. The M^/L values cal- 
culated through BdJ's colour- A/* /L relation are more scat- 
tered than those calculated from inferred SFH's, likely be- 
cause BdJ's colour- M, /I/ relations were derived for gaseous 
discs, and are not ideal for early- type dwarf galaxies. 

Immediately apparent in both plots is an offset in M, /L 
between dl and dE galaxies. The dE's have a higher M^/L 
on average than dl's, as expected from their older, redder, 
stellar populations. Median values for M^/L according to 
type are given in Table [5] 

The abnormally high Mt/ L ratio for Ursa Minor dSph 
(UMi) seen in the plot for BdJ values reflects its unusually 
red colour (see Table [l]). UMi's redness makes it an outlier 
also in Fig. [S] However, the Mt/L estimated from the SFH 
approach, which does not require colour information, yields 
a reasonable value in agreement with the observed trend for 
other Local Group galaxies. 

We compare in Fig.[2]the stellar masses computed from 
the M*/L ratios of the two methods. There is good agree- 
ment, with the exception of UMi. 

The stellar masses and stellar mass ratios computed 
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Table 1. Data for the Local Group Dwarf Galaxies 



Typo" Mt/Lv 



log M, " 



vf 



log Z" 



(S - V)„ 



M 33 


dl 


1.2 


9.55 


8.51 


0.37 ± 0.01 


135 ± 


15 


1800 


-2.7 


0.24 


0.61 


2,4,6,10,13, 
16,18,19,20 


LMC 


dl 


0.7 


9.19 


8.14 


0.17 ± 0.02 


72 ± 


7 


500 


-2.3 


0.26 


0.45 


1,7,11,16,19 


NGC 55 


dl 


0.8 


9.04 


8.15 


0.17 ± 0.05 


86 ± 


3 


1404 


- 


0.18 


0.49 


14,16 


SMC 


dl 


0.8 


8.67 


7.63 


0.15 


60 




420 


-2.9 


0.046 


0.48 


7,11,16,17,19 


M 32 


dE 


1.3 


8.66 


12.09 


- 


96 ± 


19 


2.5 


-2.8 


- 


0.93 


7,14,16 


NGC 205 


dE 


1.4 


8.65 


8.54 


-0.39 


89 ± 


15 


0.39 


-2.2 


- 


0.64 


7,14,16,19 


IC 10 


dl 


0.9 


8.30 


8.37 


-0.44 ± 0.08 


47 ± 4 


98 


-3.0 


0.06 


0.53 


7,14,15,19 


NGC 6822 


dl 


0.8 


8.23 


7.91 


-0.46 ± 0.07 


51 ± 


3 


140 


-2.9 


0.021 


0.55 


7,14,16,19 


NGC 147 


dE 


1.6 


8.17 


8.12 


-0.34 ± 0.07 


44 ± 


10 


0.005 


-2.8 


- 


0.75 


7,14,16,19 


NGC 185 


dE 


1.0 


8.14 


8.51 


-0.51 ± 0.08 


46 ± 4 


0.13 


-2.5 


- 


0.74 


7,14,16,19 


NGC 3109 


dl 


0.8 


8.13 


6.84 


0.09 


68 ± 4 


820 


-3.4 


0.02 


0.45 


7,14,16 


IC 1613 


dl 


1.0 


8.03 


7.42 


0.05 ± 0.24 


37 ± 


3 


58 


-3.1 


0.003 


0.69 


7,14,16,19 


WLM 


dl 


0.9 


7.65 


8.36 


-0.04 ± 0.03 


23 ± 


2 


63 


-3.1 


0.0011 


0.58 


7,14,16,19 


IC 5152 


dl 


0.5 


7.54 




- 


38 ± 4 


67 


-3.1 


- 


0.32 


7,14,16 


Sextans B 


dl 


0.8 


7.54 


- 


- 


38 ± 


9 


44 


-3.8 


0.002 


0.49 


7,14,16,19 


Sextans A 


dl 


0.5 


7.43 


6.82 


-0.20 ± 0.09 


33 ± 


2 


54 


-3.6 


0.002 


0.34 


7,14,16,19 


Fornax 


dE 


1.2 


7.27 


7.29 


-0.39 


20 ± 


3 


0.7 


-2.9 


- 


0.61 


7,14,16,19 


EGB 0427+63 


dl 


1.0 


6.96 


6.99 


- 


33 ± 


10 


17 


- 


0.0004 


0.55 


14,16 


And I 


dE 


1.6 


6.85 


7.01 


-0.47 ± 0.02 






0.096 


-3.1 


- 


0.75 


7,14,16,19 


Pegasus 


Tr 


1.0 


6.83 




- 


17 ± 


3 


3.4 


-3.7 


0.0003 


0.55 


7,14,16,19 


And VII 


dE 


0.9 


6.67 


7.10 


-0.54 


- 




- 


-3.2 


- 


0.51 


3,7,8,16 


Leo I 


dE 


0.9 


6.66 


7.57 


-0.85 ± 0.05 


17 ± 


3 


0.009 


-3.1 


- 


0.76 


7,14,16,19 


And II 


dE 


1.0 


6.63 


6.75 


-0.51 ± 0.02 


18 ± 


5 


- 


-3.2 


- 


0.54 


5,7,14,16 


UKS 2323-326 


dl 


0.6 


6.51 


6.50 


- 


- 




- 


- 


- 


0.39 


14,16 


GR 8 


dl 


0.7 


6.42 


7.49 


-0.81 ± 0.07 


21 ± 


6 


9.6 


-3.7 


0.0007 


0.34 


7,14,16 


Antlia 


Tr 


1.0 


6.41 


6.84 


-0.37 


12 ± 


3 


0.72 


-3.6 


0.00028 


0.52 


7,14,16,19 


Sag DIG 


dl 


0.4 


6.36 


6.44 


- 


14 ± 4 


8.6 


-4.0 


0.000067 


0.28 


7,14,16 


And III 


dE 


1.8 


6.27 


6.70 


-0.78 ± 0.04 


- 




0.09 


-3.4 


- 


0.54 


7,14,16 


Leo A 


dl 


0.5 


6.26 




- 


18 




7.6 


-3.8 


0.000032 


0.31 


7,14,16,19 


DDO 210 


Tr 


0.9 


6.25 


7.32 


- 


13 ± 


3 


2.7 


-3.6 


- 


0.25 


7,12,14,16 


And VI 


dE 


0.5 


6.17 


6.55 


-0.53 ± 0.03 


- 




- 


-3.4 


- 


0.34 


3,7,8,9,16 


Leo II 


dE 


1.6 


6.16 


7.16 


-1.05 


13 ± 


2 


0.03 


-3.3 


- 


0.63 


7,14,16,19 


Phoenix 


Tr 


1.8 


6.10 


- 


- 


17 ± 


3 


0.17 


-3.6 


- 


0.59 


7,14,16 


Seulptor 


dE 


1.7 


6.08 


7.31 


-0.77 ± 0.02 


12 ± 


2 


0.09 


-3.2 


- 


0.68 


7,14,16,19 


Tucana 


dE 


1.6 


5.97 


6.76 


-0.91 ± 0.08 






0.015 


-3.4 


- 


0.67 


7,14,16,19 


Draco 


dE 


1.8 


5.96 


7.07 


-0.99 


16 ± 


3 


0.008 


-3.7 


- 


0.92 


7,14,16,19 


Sextans 


dE 


1.6 


5.93 


6.28 


-0.51 ± 0.11 


12 ± 


2 


0.03 


-3.6 


- 


0.65 


7,14,16,19 


LGS 3 


Tr 


1.0 


5.86 


6.68 


-0.85 


17.5 ± 


0.6 


0.2 


-3.4 


- 


0.70 


7,14,16,19 


UMi 


dE 


1.9 


5.75 


6.63 


-0.79 ± 0.14 


20 ± 


2 


0.007 


-3.6 


- 


1.27 


7,14,16,19 


Carina 


dE 


1.0 


5.71 


6.38 


-0.82 


13 ± 


3 


0.001 


-3.5 


- 


0.64 


7,14,16,19 


And V 


dE 


1.1 


5.60 


6.67 


-0.84 


- 




- 


-3.6 


- 


0.57 


3,7,16 



"Morphological type simplified from iMated il99Sl) . lvan den Berghl i2000l) . and lGrebel et al.l J2003l ) as described in the text. 

Our derived values of M^/L\^ for the "b" set as described in the text. 
^Stellar mass, in units of Mq , derived as described in text. The typical uncertainty: 0.17 dex. 

Central stellar surface mass density, in units of Mq kpc~ , derived as described in text. The typical uncertainty: 0.17 dex. 
^Exponential radius in kpc. 

^Estimated velocity, in kms~ , defined as the maximum of ^rot and 2t7. 
^H I gas mass in units of 10 Mq taken from lGrebel et al.l J2003l ). 

Following fj2] log(Z/0.019) — [Fe/H]. Typical uncertainty: 0.2 dex (Eva Grebcl, private communication). 
^Current star formation rate in Mq year" . 

^De-reddencd B ~ V colour. Reddening values are from lSchlegel et al.l J1998I) except for IC 10 \vhose source is lRicher et al.l i200ll) . 

References. — (1) lAlves fc NelsonI i2000l '): (2) iBaggett et al.l i 19981 '): (3) ICaldwelll )1999D : (4) ICorbelli fc Saluccil i2000l '): (5) ICote et al] i 19991 '): (6) 
lEngargiola et al.l i2003l '): (7) iGrebel et "all i2003l '): (8) iGrebel fe Guhathakurtal )1999|) : (9) iHopp et all )1999|1 : (10) Ijacobssoni |1970D : (11) iLarsen et "ahl 
i2000l ): (12) iLee et al.l i 199911 : (13) iLee et al.l )2002ll : (14) iMated )1998|1 : (15) [Richer et ail )200lll : (16) ISchlegel et al.l )l99el ): (17) IStanimirovic et al.l 
J2003); (18) lTiede et al.l iJOoi) : (19) Ivan den Berghl )2000D : (20) IZaritskv et al.l h989D : 



through BdJ are generally slightly higher than those com- 
puted via SFH. This may be explained by the fact that BdJ 
calibrated their colour- M, /L correlation under the assump- 
tion of maximal discs, making the BdJ Mt,/ L ratios upper 
limits. 

The distribution of stellar masses, M* , against absolute 
V magnitude. My, is sho-wn in Fig. [S] The solid line has a 
forced slope of -0.4, while the zero-point is adjusted to fit 
the data. Once again, the M* values for the dl galaxies tend 
to fall belo-w the dE's as a result of the younger populations 
in dl's. 

Both methods of calculating M./L yield comparable 
values of log M, within 3% (see Table [2}. Thus, the ob- 



served luminosity scaling relations (e.g. E — L, [Fe/HJ-L, 
Re-L, [0/HJ-L, V-L, a — L) will remain intact when trans- 
lating them to physical scaling relations. 



2.1.3 Choice of Methods 

Given the incompleteness of our original data bases, we es- 
tablish various guidelines in Table |3] to assign final M, /L 
values per galaxy. 

We can use a median M^/L value for all the galax- 
ies according to Hubble type, or assign the M^/L for each 
galaxy based on the BdJ or SFH methods. The first option 
recognises the uncertainties of both methods of calculating 
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Figure 1. The BdJ and SFH stellar niass-to-light ratios, in the V-band, as a function of My- 



M,/L (namely, the uncertainties in the SFH, and the un- 
certainties in extrapolating BdJ's colour- M./L correlation 
down to dwarf galaxies), while the second method uses the 
full available data to make educated guesses for M^,/ L. We 
have thus created two data sets of stellar data according to 
these methods, calling the first "a," and the second "b." 

For the "a" set, we assign the average of the SFR and 
BdJ median values of the M»/L from Tabled For the "b" 
set, we first shift the BdJ M^/L ratios down so that their 
median values per Hubble type match those from SFH. We 
then assign their M*/L ratios according to our prescription 
in Table El 

The "a" and "b" sets are consistent with each other 
within their errors, and our choice of "a" or "b" makes no 
difference to our conclusions. In the subsequent analysis, we 
employ the "b" set since it utilises all the available data to 
calculate M^/L, and preserves any intrinsic scatter. 
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Figure 2. Difference in stellar masses computed via BdJ and 
SFH. 



Table 2. 


Median M, 


/L 


for different 


galaxy types 












S-B 


ind 


- 


V-B 


md 




SFH 


BdJ 


SFH 


BdJ 


dl 


0.5 


0.6 




0.7 


0.8 


Tr 


0.7 


0.8 




1.0 


1.0 


dB 


1.5 


1.1 




1.6 


1.3 



2.2 Surface Density 

The luminosity profiles of the dl and dE galaxies have been 
fitted (M98 and references therein) with an exponential 
function whose inner extrapolation provides the central sur- 
face brightness So, which we translate to a stellar surface 
density, /x, (M0kpc~^) using the M*/L derived above. 

M 32 has a very high central surface brightness of 11.5 
magarcsec"^ (GGH03), yielding log /x, = 12.1 (M© kpc"^). 
This is four orders of magnitude higher than other galax- 
ies of similar luminosity. M 32 ind eed belongs to the rare 
class of compact d warf galaxies fe.g. lZiegler fc Benderlll998l : 
iMieske et al.ll2005l . and references therein) . We thus exclude 
M 32 from analyses involving ^«. However, we retain it for 
other tests. 
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Figure 3. Comparison of stellar masses with the absolute magnitude My . The BdJ and SFH computed values are on the left and right 
respectively. The solid line has a forced slope of -0.4 and a free vertical offset. 



2.3 Velocity 

The velocity data for the dwarf galaxies consist of a ro- 
tational velocity (corrected for inclination), K-ot, and/or a 
projected central velocity dispersion, a. The dl's typically 
have a rotation velocity that is larger than the velocity dis- 
persion (exceptions ar e Leo A, GR 8 and Sag DIG - see 
ICourteau et al.ll2007bl ). whereas dE's only have a measur- 
able velocity dispersion. 

We wish to define as a global parameter a character- 
istic velocity V that will represent in both cases the depth 
of the potential well, which for dwarf galaxies is believed 
to be dominated by the dark halo. The measured rotation 
velocity in the disky dl's is taken as is, assuming that it 
approximates the virial halo velocity and is a near lower 
bound to the maximum rotation velocity. We then ask, what 
does a dE's velocity dispersion tell us about its potential 
well? What would have been its rotation velocity if it had 
been disky in the same halo with the same stellar mass? 
To answer these questions, given the tightness of the TuUy- 
Fisher (TF) relation, we assume that the halo mass scales 
with the stellar mass in the same way for both dE's and 
dl's (refer to i]4.1|l and adopt the following algorithm: take 
V = max{ V^ot, Xcr} (or the appropriate value when only one 
of Vrot and a is available) where X is chosen to minimise the 
scatter in the TF relation (log V vs. log Af«). We find that 
X = 1.92 whi ch is consistent with th e V-a relation for large 
disc galaxies IjCourteau et al.ll2007a ). but is larger than the 
X = \/2 solution for an isothermal sphere. Had we chosen to 
minimise the baryonic TuUy-Fisher (BTF) relation instead, 
we would have found a lower value of X but slightly larger 
scatter than for the TF relation for these dwarf galaxies (see 
also i|4.1l for further motivation to minimize the TF instead 
of the BTF relation). The choice of the isothermal solution 
X — \/2 would yield a TF relation that is slightly steeper 
(in log V vs. log M,), but the galaxy types still form a sin- 



gle relation. Furthermore, the results of our PCA analyses 
(described in 331) are statistically identical for both X = \/2 
and X = 1.92. 

In summary, our algorithm sets V = Xa for the dE 
galaxies, and V — Vrot for the dl galaxies (except for GR 8, 
Leo A and Sag DIG where a is greater than the maximum 
Kot of the ISM). 



2.4 Radius, Baryon Mass, Metallicity and SFR 

We use the exponential scale length scaled with the distance 
estimates provided in M98 as the characteristic stellar radius 
of galaxies. We also use the H I mass provided by GGH03 
and our derived stellar mass to estimate a baryonic meiss. 
For most of the dE's the quoted H I mass is an upper limit. 
Iron abundances [Fe/H] are taken from GGH03, which 
are based on spectroscopic and photometric measurements 
of red giant branch stars. GGH03 use the mean abundances 
of the old stellar populations which can be consistently mea- 
sured across the galaxy types. Oxygen abundances are avail- 
able only for galaxies with current star formation, so besides 
the dl's, only four dE's and one Tr galaxy have available oxy- 
gen data. Thus, we use the iron abundances to trace the to- 
tal metallicty Z and assume a linear relation between them. 
Recalling that [Fe/H] is logarithmic, we then define Z as: 



log{Z/ZQ) = [Fe/H] 
where we ad opt Za 



(1) 



iGirardi et"aLll2000l : ICarraro et al 



0.019 (lAnders fc Grevessi Il989l : 



20011 ). 



The current star formation rate (SFR or M,) data for 
the dl galaxies, estimated from observed Ha flux, are taken 
directly from M98 or vdBOO. dE galaxies have very little or 
no current star formation. 
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Table 3. M,/L assignments for data set "b" 



Case 1 



dl average M^ / L ratios from BdJ and SFH 

dE M,/L from SFH 
Tr M,/L from SFH 



\vhichcver datum exists 
wiiiclrcver datum exists 



AI^/L from data set 
AI^/L from data set 



Note. — Tire eases are: 1. enough data exist to ealculate AI^/L from botli SFH and BdJ methods; 2. 
data exist for only one method of computing A/*/L; 3. not enough data exist for either SFH or BdJ 
methods. (For Tr galaxies, only Case 1 applies.) 



2.5 Error bars 

Due to the uncertainties in the determination of M,/ L and 
the heterogeneous data base, error bars are not taken into 
account in the subsequent analysis. However the typical er- 
ror bar size relative to the scatter in distribution of the 
galaxies is shown in the upper left corner of all figures. The 
estimated error bars for logM* and log/x, reflect the dif- 
ference between the SFH and BdJ calculations Mt,/L, and 
are on average 4% and 7% of their respective ranges. The 
error bars are typically 4% for log V and 10% for log Z (Eva 
Grebel, private communication). No estimates are provided 
for the uncertainty in SFR. 



3 METHODS OF ANALYSIS 

Throughout this analysis, we study the relations between 
the dwarf galaxy structural and SF parameters using linear 
regression and principal component analysis (PC A). 

For the two-dimensional analyses, we investigate the 
correlation in log space between pairs of the structural and 
SF quantities. We plot the distribution of each pair, obtain 
the best-fitting line by the method of bisector least-squares, 
and calculate the st andard Pear s on lin ear correlation coeffi- 
cient r according to lPress et al.l (|l992h . 

We use PCA as a tool to quantify the shape of the distri- 
bution of the galaxies in parameter space. Given a data set 
with m parameters (such as stellar mass, metallicity, etc), 
the PCA quantifies the distribution of the data along orthog- 
onal basis vectors of the parameter space. PCA outputs the 
vectors V^, fc = 1, m and their eigenvalues Dk. Appendix B 
describes the derivation of Vfe and Dk (or Dk,k as in the 
Appendix) . 

The following is an interpretation of the eigenvectors 
and eigenvalues: 

The eigenvectors, which are orthonormal, lie in the di- 
rections of greatest variance in the data. The ratios of the 
eigenvalues are a measure of the relative strength of the vari- 
ance along each corresponding eigenvector. For example, if 
the distribution of the data were a line in a parameter space 
of 3 parameters, the vector Vi with the largest eigenvalue 
D\ will lie along that line, while the other two vectors will 
be orthogonal to it. Furthermore, the eigenvalue D\ will be 
much larger than the other two eigenvalues, which will be 
similar in size, reflecting the extended distribution of the 
data along Vi, and relatively little extension in the other 
two directions. Such a distribution is characterised by one 
primary parameter. If the distribution of the data is a plane 
in 3-D space, two of the eigenvectors will lie in the plane 
and will have comparably large eigenvalues, while the third 
vector, lying orthogonal to the planar distribution will have 



a relatively small eigenvalue. Such a distribution is charac- 
terised by two primary parameters. 

In fact, for 3-D space, we can more generally associate a 
"linear" distribution with "prolate" (Di:D2 3> D2:Dz), and 
a "planar" distribution with "oblate" {Di:D2 ^ D2'-D3). 

In order to meaningfully quantify the shape of the data 
distribution, we standardise the data before performing the 
PCA. Standardisation divides the data set by the standard 
deviation of the parameters, so that the eigenvalues are in- 
dependent of parameter range. In other words, standardis- 
ation gives us a shape for the distribution with the ranges 
of all the axes on equal footing. Without standardisation 
the shape of the distribution will change from being linear 
to planar, depending on their relative scalings. In fact, the 
scaling relations themselves can be derived directly from the 
elements of the first eigenvector of the unstandardised anal- 
ysis. On the other hand, standardisation will output eigen- 
values that are independent of the scalings, and therefore 
insensitive to whether we choose for instance to use V^ in- 
stead of V. Therefore, throughout this study, we present 
the results of standardised PCA along with the standard 
deviations of each parameter. (The scaling relations may be 
recovered by dividing the components of the first principal 
eigenvector Vi by the corresponding standard deviations.) 

We caution the reader that standardisation assumes a 
complete, unbiased representation of all galaxy types for all 
considered parameters. Given the depth of current imaging 
surveys, we suspect that only much fainter LG dwarfs are 
truly missing, and our assumption of standardisation is rea- 
sonable. Performing our analysis using strictly unstandard- 
ised PCA would yield tighter fundamental lines in every case 
due to the large range of M, relative to the other parame- 
ters. Standardisation is nonetheless preferred for the reasons 
stated above. 

In the following sections, we present our 2D and PCA 
analyses first for the structural parameters, and then for the 
SF parameters. 



4 RESULTS: THE STRUCTURAL 
PARAMETERS 

4.1 2D: Velocity and Mass 

Fig- m shows the correlation of characteristic velocity versus 
A/* for LG dwarf galaxies with V oc M°-^'^^°-°^ and r = 
0.94. Recall that V is constructed from Kot and a in such 
a way as to minimise the scatter in this scaling relation. 
So naturally, the dl and dE galaxy types seem to lie on 
the same relation with uniform spread. However, had we 
used V — (\/2 or \/3)o" instead of ~ 2a the dl's and dE's 
would still appear to lie on the same relation. Note also that 
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Figure 4. Left: Circular velocity vs. stellar mass for LG dwarf galaxies. Right: Circular velocity vs. "baryonic mass" (Af^ar = Af»+AfHi). 
We show bisector fits for all the galaxies (solid), and separately for the dE (dotted) and dl (dashed) types. The dash-dotted line delineates 

each figure. 



ilaxies below M« = 10^'^Mq. The Pearson correlation coefficient r is shown at the top and the bisector fit slope s at the bottom of 



the smallest dwarfs cluster just above V ~ lOkm s~ and 
could even be described b y a constant V (dash-dotted line) 
which lDekel fc Wool l|2003h explain might be due to radiative 
feedback. 

The LG dwarf galaxies fall below K ;S 90kms"^ The 

V — M» r elation of LG dl's is consistent with that of 
iMcGaughl (|2005l '). with a slope of 0.26 ± 0.02 or M, oc 
y3.8±o.3^ which is not significantly different from M» oc V"^ . 

Adding the gaseous mass from GGH03 {Mi,ar = Mt + 
Mm) yields a baryonic TuUy-Fisher (BTF) slope for the LG 
dl's that is si gnificantly shallo wer (M, oc V^3.3±o.2-j ^^ig,^ that 
reported by JMcGaughlbOOSl). but cons i stent with the BTF 
relations of both BdJ and iGeha et~aLl (|2006l ) (Fig. H left). 
We also find that for dFs with log V > 1.5, we find that the 
scatter in the LG BTF relation is reduced significantly com- 
pared to Fig. |4] (right). The se dl's lie almos t precisely along 
the BTF re lation foun c l by l|Mc"Gaughll2005l 'l . with compara- 
ble scatter. iMcGauehl (J2005|) describes a break in the distri- 
bution of circular velocity versus stellar mass (or total light) 
for galaxies with Vc ^ 90kms~^. This break however disap- 
pears if Vc is plotted against total (luminous + gaseous) 
mass. 

The shallower overall slope of the dl BTF relation is 
likely due t o the d ifferent prescriptions for calculating M,/L. 
iMcGaughl (|200a ) used an optimal empirical relation (see 
their section 3.3), while we used stellar population mod- 
els, as did BdJ (see ^2TT |> . IMcGaughl l|2005h showed that this 
difference in methods is enough to account for the shallower 
BTF slope found by BdJ (which is consistent with ours). 

The contribution of the H I mass to dE galaxies is 
negligible and we find two distinctly offset BTF relations 
for dl and dE galaxies. (This is true even had we used 

V — \/2cr for the dE's.) This off set would have been larger 
had we used the suggestion of IPfenniger fc Revaa (|2005l ) 



that Albar = Af, + 3Mhi. This suggests that although the 
BTF may more fundamentally describe the kinematics of 
dl's, the stellar relation (i.e. the TF) may be a more funda- 
mental description of the entire dwarf population (since the 
same TF describes all galaxy types). (Recall also that had 
we minimised the scatter of the BTF in determining the X of 
i]2.3l the scatter of the BTF would have been slightly larger 
than the minimised scatter of the TF.) Using baryon mass 
instead of stellar mass also increases the offset between the 
types in the i?* — Afbar, pt. — Mhm, and Z — Mbar relations. 



4.2 2D: Radius and Mass 

Fig. [S] shows the correlation between exponential radius and 
stellar mass fitted as R, oc M°-^^^°-°'^ , with r = 0.85. The 
galaxy types populate different regimes of the relation, with 
dl galaxies being bigger and more massive than dE's, yet 
they are both described by the same R* — M* relation. This 
is interesting in itself as it may indicate a similar physical 
origin to the relation, (such as halo mass), which may simul- 
taneously be responsible for their differences. 

Shon ct al. (2003) have plotted the median of the dis- 
tribution of the sizes of the SDSS galaxies as a function of 
stellar mass (see their Fig. 11). The sizes are measured for 
the z-band and are given as Sersic half-light radii. If we ne- 
glect bandpass differences, we can make a rough comparison 
of the LG data with SDSS data. Gonsidering that the expo- 
nential half-light radius is R1/2 = 1.687?*, we find that the 
LG dE's appear to be the low-mass extensions of the early- 
type distributions of sizes for the SDSS galaxies. The SDSS 
early-type galaxies follow the relation 7?« oc M*'^^ while 
LG dE's have a shallower R-^ oc M°-^'^^°-°^ . This shallower 
slope is consistent with the lower end of the size-luminosity 
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Figure 5. Exponential radius vs. stellar mass. We plot the two- 
dimensional linear regression on all the galaxies (solid), and sep- 
arate fits to the dE (dotted) and dl (dashed) types. 



relation of iGraham fc WorlevI (|2008l ) while we extend their 
relation to lower masses. 

However, while the LG dl's overlap the /?* — M* rela- 
tion of the SDSS late-type galaxies, they extend the relation 
which a steeper slope of 0.41 ± 0.04 compared to the slope 
of 0.15 at the low- mass end of the SDSS sample. 

We can derive a stellar size-mass, i?* — M*, relation for 
large spiral galaxies from the comp ilation of scaling param - 
eters for 1300 local disc galaxies bv lCourteau et al.l (l2007al ). 
We use the /-band exponential scale radii and luminosi- 
ties, V — I colours , and colour- Aft/L relations of BdJ (see 
iDutton et al.l 120071 ) to derive a Rt — M« relation for large 
galaxies. We find that the dl's overlap with the R, — M, 
relation of bigger disk galaxies and with comparable slope. 

4.3 2D: Surface Mass Density and Mass 

The dependence of radius and stellar mass is also expressed 
via the relation for stellar surface density {/it, ex Mt./R'l) 
and stellar mass. Fig. |6] shows a correlation between /i* and 
M* that spans 4.5 decades in M, and 3 decades in /i,. The 
best-fitting relation is /x. oc M°-^'^^°-°^ with r = 0.81. 

With the exception of WLM and GR 8, dl galaxies have 
a lower mean surface mass density than dE's. Since both M* 
and fi, depend on M^/L in the same way, plotting /i, against 
M* or L, as previously done by GGH03, yields identical 
offsets between the dE and dl types. These best fit lines 
have the same slope of 0.74 ± 0.05, but are offset by about 
0.5 dex from the centre of the distribution. 

iKauffmann et al.l (|2003bl . see their Fig. 7) showed an 
analogous distribution in the plane of surface density (within 
the half-light radius) versus stellar mass for 123,000 SDSS 
galaxies with A/* > lO^M©. They find a weak systematic 
dependence of about /x, oc Af°'^ (our eye-ball estimate) for 
the HSBs and ^* oc M°-^^ (their quote) in the LSB regime 



Figure 6. Central stellar surface density vs. stellar mass. We 
plot the two-dimensional linear regression line of all the dwarf 
galaxies (solid), and separate regressions for the dE (dotted) and 
dl (dashed) types. 



with 10* < AT. < 3 X 10^°Mq. After considering the /x. - 
^e conversion, the LG dwarf galaxies extend this relation 
down to 8 X 10^ M© and together, their /x* — A/, relation is 
con sistent with /j« oc M°'^ . 

IShen et al.l l|2003h showed the median of the same distri- 
bution, but separated into late and early-type distributions. 
They show that the early-types tend to have higher surface 
density, consistent with our finding for LG galaxies. In fact, 
the LG dE's and dFs again seem to extend early and late- 
type surface density relations to lower masses. (Our eye-ball 
estimate of their slopes: ~ 0.65 for both late- and early- 
typ es, compared to s — 0.6 1 ± 0.05 for LG dwarfs). 

IVaduvescu et al.l (120051 ) studied a sample of 34 dl galax- 
ies and provided central surface brightnesses and K and J- 
band magnitudes for these galaxies using a hyperbolic secant 
(sech) function to fit their surf ace brightness p r ofiles. Us- 
ing the colour- Af,/L relations of iBeU fc de Jond (|200ll ). we 
converted their sech magnitudes and central surface bright- 
nesses to stellar quantities, and compared their Af«-/x« re- 
lation to that of the LG dwarf galaxies. We find that the 
Mt-fj., relation is significantly tighter (r = 0.86) than their 
magnitude relation (Afx,sech-Mo,sech, r = 0.71). Secondly, we 
find that their dFs lie neatly along the same Mt-jit relation 
as the LG dwarf galaxies, even though the surface brightness 
profiles of the LG dwar fs were fit to exponentia l functions. 
We plot the data from IVaduvescu et al.l (|2005l ) , which we 
converted to stellar quantities, together with the LG data 
in Fig. [7] 



4.4 PCA: The Fundamental Line in Structural 
Space 

We use PCA as outlined in |3] to quantify the shape of the 
distribution of the LG dwarf galaxies in the (log) space of the 
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Figure 7. The Mt-fi* r elation of the LG dwarf galaxies plotted 
together with the dl's of lVaduvescu et alj 1120051 ) (open circles). 



structural parameters (M«, R,, V). We present the results 
graphically in Fig. |8] and quantitatively in Table |4] 

Fig. [S] shows the axes Yj (defined in Appendix B) as 
the coordinates of the data vectors projected onto Vj, for 
j = 1...3. Each panel represents different views along the 
three dimensions of the distribution. On visual inspection, 
the shape of the distribution is linear in Fig. [8] 

We list the components of the three Vj vectors in Table 
m The first vector points in the direction of the fundamental 
line in this scaled basis. The other two vectors point in the 
directions orthogonal to the line. By construction, the com- 
ponents of the first principal vector Vi are ~ l/\Ai where 
n is the number of dimensions (in this case 3) since stan- 
dardisation places all the data along a diagonal. The vector 
components are in the basis of the original four parameters 
scaled by their standard deviations, from which the scaling 
relations are derived. (In the unstandardised analysis, the 
scaling relations are recovered directly from the components 
of the first eigenvector.) For example, if the first eigenvec- 
tor is Vi, and its j'th component is Vi(j), then the scaling 
relation between the parameters log M* and log V is 



(2) 



Vi{M^)/aM, _ 0.60/1.22 
yi{V)/av ^ 0.58/0.32 

yielding log V oc (0.27±0.01) log M, which is consistent with 
the 2-D fit in i|4.1l The scaling relation between log M* 
and log Rt is listed in the note under Table |4l and is also 
consistent with the 2-D fit in ^4.21 

The eigenvalues associated with the principal vectors 
quantify the shape of the distribution of the data. The eigen- 
value Eissociated with the first principal vector is 90.1% of 
the sum of all the eigenvalues. The strength of the second 
eigenvalue is only 8.4%. The variance along the first princi- 
pal vector is more than 10 times greater than the variance 
along the second principal vector. The ratio of the second-to- 
third eigenvalues is only 5.6. In other words, the extensions 



Figure 8. The fundamental line in structural space. The axes 
lie in the directions of the Vj eigenvectors defined in the text 
and the data are projected onto them. The components of the 
eigenvectors are given in Table |4] 



of the data along the second and third principal vectors are 
more similar to each other than they are to the extension 
of the data along the first principal vector. This describes a 
linear or prolate distribution, according to the criteria in ^ 
and is characterised by one primary parameter. 

The standard method of displaying fundamental lines or 
planes (edge-on) from PGA is to plot one parameter against 
a linear combination of the other parameters. So the LG 
data would be plotted as log A/» versus 

(2.92 ± 0.31) log V + (0.72 ± 0.32) log R, + (3.35 ± 0.23) (3) 

from the third principal vector. This is equivalent to display- 
ing only the upper-left panel in Fig. [S] but rotated so that 
the log M* lies along the y-axis, and scaled so that the data 
lie on the y = x line. Such a presentation of the data, with 
its confusing mix of parameters on its x-axis, is not spatially 
intuitive as it fails to display their distribution from other 
angles. We feel that our method of displaying the output of 
PGA is more useful since the data are displayed along all 
the principal axes of their distribution, giving the viewer a 
more intuitive feel for its shape. 

We also performed PGA on an "alternate" structural 
parameter space where we have replaced 7?« with fj,^ in the 
above structural space. We summarise the results in Fig. [5] 
and Table [5] The distribution of the galaxies is triaxial in 
this space since the ratios of the eigenvalues (listed in the 
table) are comparable. However, since the first eigenvalue 
is at least (8.0) times greater than the other eigenvalues, 
meaning that the galaxies are at least 8 times more extended 
along one axis than the along the other axes, we say that 
the galaxies lie in a linear distribution, characterised by one 
primary parameter, in this alternate structural space. 
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Table 4. PCA results in structural space 



log A/* 



logV 



logfl. 



Eigcn- Ratios 

values (%) 



± 0.01 0.58 ± 0.01 0.56 ± 0.01 90.1 ± 3.2 



V2 -0.23 ± 0.08 -0.54 ± 0.07 0.81 ± 0.03 
V3 0.77 ± 0.02 -0.61 ± 0.06 -0.19 ± 0.08 
ax 1.22 0.32 0.41 



1.5 ± 0.5 



}10.7± 3.9 
}5.6± 2.7 



Note. — The scaling relations from the principal vector: 
logV ex (0.27 ± 0.01)logM,, log H, cc (0.36 ± 0.01)logM, . 
Standard projection: 
logM, = (2.92 ± 0.31)log V+{0.72 ± 0.32)log fl.+(3.35 ± 0.23). 
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Figure 9. Projections of the LG dwarf data in an alternate struc- 
tural space. The components of the PCA eigenvectors are given 
in Table [5] 



5 RESULTS: THE STAR FORMATION (SF) 
PARAMETERS 

5.1 2D: Metallicity and Mass 

Fig. 1 101 shows metallicity versus Af, . The best-fitting global 
scaling relation is Z x M°*°*°°* with r = 0.76. The higher 
metallicity of dE galaxies relative to the dl galaxies of the 



Table 5. PCA results in alternate structural space 



log M, 



logV 



Eigen- Ratios 

values (%) 



Vi 0.60 ± 0.01 0.54 ± 0.02 0.58 ± 0.01 87.5 ± 3.8 1 „ „ , „ 7 

V2 0.24 ± 0.06 -0.82 ± 0.02 0.51 ± 0.05 11.0 ± 3.7 |_' , ^' 

V3 0.76 ± 0.02 -0.17 ± 0.07 -0.63 ± 0.04 1.6 ± 0.5 -'' 

ax 1.19 0.71 0.32 

Note. — The scaling relations from the principal vector: 
log;/, X (0.66 ± 0.03)logM, , log V oc (0.27 ± 0.01)logM, . 
Standard projection: 
logM, = (0.35 ± 0.15)logM.-l-(3.00 ± 0.24)log V+(0.28 ± 0.37). 



same luminosity has been known for some time (see for ex- 
ample M98) , and this offset persists even for old stellar pop- 
ulations (GGH03). Not surprisingly, we see the same phe- 
nomenon when using M« instead of luminosity. Their respec- 
tive correlations are much stronger when fitted separately 
(dE's: 0.93; dl's: 0.88)^ 

For SDSS data. lKauffmann et all ()2003bl ') reported Z ex 
j^o.4-0.5 (depending on the metallicity tracer) for Af, < 
3 X lO^^Mp), at the brig ht end of the dwarf galaxy regime 
(JKauffmann et al.ll2003bf ). This is consistent with our global 
and type-dependent slopes. 

iTremonti et al.1 (1200J) showed more specifically that 
SDSS star-forming galaxies follow the relation 12 -I- log 
[0/H] ~ 0.35 log M, for M. < 3 x IO^^Mq, with an over- 
all saturation beyond this transition mass (see their Fig. 6). 
This translates to Z oc MS'^^ at the low-mass end assuming 
that the oxygen abundances trace Z linearly. Using oxygen 
abundance data from M98, we find that the LG galaxies fol- 
low 12 -I- log [O/H] ~ 0.39 ± 0.04 log M*, consistent with 
the slope of low mass SDSS galaxies (Fig. Ill|l . Our slope is 
also a slightly better match to the low -mass end of the SDSS 
galaxies than that of iLee et alj (|2006l ) . 

However we find that the LG dwarf galaxies have oxy- 
gen abundances that are on average 0.3 dex lower than 
those of low-mass SDSS galaxies. While oxygen abundances 
of the LG dwarf galaxies are measured via spectroscopy of 
H II regions (M98), abundances for SDSS galaxies were es- 
timated using empirically and theoretically calibrated re- 
latio ns between met allicity and strong optical emission line 
flux. ITremonti et al.i (2004) notes that aperture bias can lead 
to overestimates i n abundances by about 0.1 dex. They also 
cite evidence from lKennicutt et al.l l|2003l ) that "strong line" 
methods such as theirs may overestimate the true abundance 
by as much as a factor of 2. Taking these overestimates into 
account, the LG dwarf galaxies coincide with the lower end 
of the SDSS distribution and extend the mass range down 
by about two magnitudes. 



5.2 Star Formation Rate and Mass 

Fig.[T2lshows the SFR versus A/» for dFs only (and two Tr's; 
dE's have very little or no star formation). The best-fitting 
scal ing relation is M, o c Mj-^ *^^""^ with r = 0.96. 

iBrinchmann et al.l (I2OOJ) (see their figure 17) showed 
that the low-mass SDSS galaxies follow a SFR-stellar mass 
relation with a log slope of ~0.7, noticeably shallower than 
the LG dwarf galaxies. Moreover, the LG dwarf galaxies 
fall significantly below (almost one order of magnitude at 
least) the SDSS distribution. Adjusting for the IMF differ- 
enc es in the SFR det e rmina tions of M98 (Salpeter IMF) 
and IBrinchmann et al.1 (|2004J) (Kroupa IMF) would only in- 
crease the vertical offset between them. Since dwarf galaxies 
have very low metallicities, which normally means very low 
dust content, the effects of dust extinction would also be 
stronger at higher mass. Correcting for this would further 
steepen the LG relation. Underestimation of the Balmer ab- 
sorption in the Ha maps which are used to estimate SFR is 
also more likely at higher mass (Brinchmann 2005, private 
communication), and correcting for this would also lead to 
a steepening of the relation. The different slope and vertical 
offset between the SDSS and LG galaxies seems to be real. 
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Figure 10. Metallicity vs. stellar mass. We plot the two- 
dimensional linear regression over all the dwarf galaxies (solid), 
and separate fits for dE (dotted), and dl (dashed) galaxies. 



Figure 12. Star formation rate vs. stellar mass. Shown is the 
two-dimensional linear regression with its slope s indicated as 
well as the correlation coefficient r. 
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Figure 11. Oxygen abundances of the LG dwarf galaxies vs. 
stellar mass. 



and present a new theoretical challenge. We discuss this fur- 
ther in 301 



5.3 PCA: The Fundamental Line of the 
Structural + SF parameters 

The above linear correlations motivate an investigation of 
the fundamental line in higher dimensional spaces than 
the structural spaces that we have already analysed. Thus, 



adding the metallicity to the structural spaces, we per- 
formed PCA on the 4-D (log) space of M* , i?* , V, Z and on 
M, , ij.,,V, Z space, and find that the fundamental line re- 
mains linear in these spaces in that the first eigenvalue is 
much larger than the other three. The line is is illustrated 
in Fig. [T3] and the eigenvalues are listed in Tables [S] and [T] 

However, a glance at Fig. [TH] reminds us that the dE's 
and dl's are offset in metallicity and suggests that they sepa- 
rately lie in linear distributions. Thus we also perform PCA 
on the dE's and dl's separately. 

For the dE's, we find that their fundamental line lies in 
the same space as the fundamental line of the whole group 
but the dE fundamental line is much tighter as indicated 
by their eigenvalues (listed in Tables [6] and [7]) , especially in 
Mt , fit,V, Z space (displayed in Fig. [T3| . 

From looking at the ratios of the eigenvalues, the dl's 
appear oblate in M,,Rt,V, Z space, and "quadraxial" in 
Af» , iJ,,,V, Z space. However noting that the first eigenvalue 
is much larger than the others in both these spaces, the dis- 
tribution of the data extends mostly in one direction, and is 
to first order linear. Hence this distribution is also charac- 
terised by one primary parameter. Recalling the tight cor- 
relation between M* and M* in i^5.2l for dFs, we add Af* 
to the parameter space and perform the PCA. We find that 
the dl's are linearly distributed in this parameter space (see 
Fig. I16p . as indicated by their eigenvalues in Table [S] 

Thus we find that, when the stellar evolution parame- 
ters Z and Mt are considered together with the structural 
parameters, the distribution of all the LG dwarf galaxies 
remains linear. In particular, the dE's and dl's separately 
form tighter linear distributions in the parameter space of 
Mt ,fJ-*,V, Z, and also of A/. , /x. , V, Z, M, for the dFs. 
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Figure 13. The fundamental line in structural space + metallicity (M«,/it, V, Z). The eigenvalues of the PCA are given in TableJT] 



Table 6. Eigenvalues for PCA in the parameter space of M« 
V, Z 



^ following summarises how we achieved those goals and our 
results. 



All galaxies dE's only dl's only 

Values (%) Ratios Values {%) Ratios Values {%) Ratios 



81.4 ± 3.9 



} 6.2± 1.6 
} 2.8± 1.1 
} 6.4± 2.5 



84.8 ± 6.4 

11.1 ± 4.5 

3.4 ± 2.7 

0.7 ± 0.4 



} 7.6± 3.2 
} 3.3± 2.9 

} 4.7± 4.4 



80.9 ± 8.9 

11.1 ± 6.0 

7.7 ± 3.7 

0.3 ± 0.2 



} 7.3± 4.0 
} 1.4± 1.0 
} 23.4±15.S 



Table 7. Eigenvalues for PCA in the parameter space of M*, ^* 

y, Z 



All galaxies 
Values (%) Ratios 


dE' 
Values {%) 


3 only 

Ratios 


dl's 
Values {%) 


only 

Ratios 


^^■^^^■^ }8.6±2.9 
1:2 ± I, > '-'^ -5 


92.4 ± 6.2 
4.1 ± 4.3 
2.6 ± 1.9 
0.9 ± 0.4 


} 22.4±23.3 
} 1.6± 2.0 
} 3.0± 2.6 


81.9 ± 8.1 

14.4 ± 7.7 

3.2 ± 1.4 

0.5 ± 0.2 


} 5.7± 3.1 
} 4.5± 3.2 
} 6.4± 4.0 



6 DISCUSSION 

We have studied the scaling relations of the dwarf galax- 
ies in physical quantities and quantified the "fundamental 
line" (FL) in the parameter space of these quantities. The 



Table 8. Eigenvalues for PCA on in the 
parameter space of Af*, fit, V, Z, M* 



dl's 


only 


Values (%) 


Ratios 


83.9 ± 7.2 






} 6.9± 3.7 


12.2 ± 6.6 




} 4.7± 3.4 


2.6 ± 1.3 




} 2.5± 1.6 


1.1 ± 0.5 




} 4.1± 2.5 


0.3 ± 0.1 



6.1 The Data 

To derive the physical stellar mass, we used a combi- 
nation of two methods to derive stellar mass-to-light ra- 
tios. These methods involved the col our- M,/L rel ations of 
iBell fc de Jong (2001) (updated bv B eU et al.ll2003l ). and SF 
histories of LG galaxies coupled to population models of 
BC03. 

We divided the six parameters (A/,, R,, jj,,, V, Z, Af,) 
into two categories. The first category is the "structural" 
quantities (Af, , /x,, fj, , V), and the second is the "star for- 
mation" (SF) quantities {Z, Af,). We treat the stellar mass 
as the independent parameter and determine the scaling re- 
lations of the other quantities with respect to it. We list 
these scalings in Table |9l Although we have used only the 
bisector fits, for comparison with other studies, we include 
in this table the slopes for the 1-D forward least squares fits. 



6.2 The Structural Scaling Relations 

We find that the scaling relations for LG dwarf galaxies ex- 
tend the corresponding relations found in larger galaxies. 
Among the structural scaling relations, the Af,-/i, relation 
of the LG dwarfs match that of the lower-mass regime in 
the SDSS, and extend the range of the relation to about 
lO^Mp). This relation i s also consistent with that found by 
IVaduvescu et alj (|2005n after converting their photometric 
data to stellar quantities. The Af,-/?, for the LG dl galaxies 
also seem to extend the correspo nding relation for large disc 
galaxies ijCourteau et alJ 2007al'). and th at of the late-type 
galaxies in the SDSS (jShen et al.l |2003|) but with steeper 
slope. 

The slope of the TF relation for the LG dwarf galaxies 
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Figure 14. Projections of the dE data in structural space + 
metallicity {Mf,fit,V,Z). The eigenvalues of the PC A are given 
in Table [3 



Figure 16. Projections of the LG dl's in parameter (log) space of 
Mt, fj,t,V,Z and M». The components of the PCA eigenvectors 
are given in Table [8] 



>- 



> 



>- 



H ' I ' I ' I 

=1 I I I I I I ^ ^ 



-o 

-2 



2 



-2 

-4 



-A 



-A^-A^Z^ 



.1.1.1 



Standardised 
M., fj,,, V, Z 



I ' I ' I 



1 1 h 1 1 1 ^ 1 1 1 1 1 1 ' 4 



■. I . I . I I- 



■' I ' I ' I '■ 



I.I.I 



Y, 



Y. 



Y, 



-4-2 2 44-2 2 4-4-2 2 4 



Figure 15. Projections of the dl data in structural space + metal- 
licity {M*, fi,t,V, Z). The eigenvalues of the PCA are given in 
Table E] 



is consistent with that found by iMcGaughl ((20051) though 
we do not find their reported break at 90 kms~^. The LG 
dwarf galaxies follow a smooth TF relation with no indi- 
cation of a mass or velocity scale. When adding the H I 
m ass, the LG dl's follow a BTF relation that matches that 
of lMcGaugiil lJ2005l ) above F = 32 km s"\ with reduced scat- 
ter. However, including all the d l's yields a shallo wer slope 
compared to the optimal slope of lMcGaughl l|2005h . but con- 



sistent with BdJ and iGeha et al.l (|2006f l. iMcGauehl (|2005h 
contended that the BTF slope is sensitive to the method 
used to determine M*/L, so our shallower slope is expected. 



6.3 The SF Scaling Relations 

The K'h-Z relation also appears to be a low-mass extension 
of the relation for the SDSS galaxies. 

The slope of the Mt-Mt relation for the LG dl's 
(1.18 ± 0.08) is significantl y steeper than that of th e SDSS 
star-forming galaxies (0.7, iBrinchmann et al.ll2004 ). Addi- 
tionally, the LG dl's fall at least one order of magnitude 
below the SDSS distribution. Correcting for IMF differences 
and observational effects such as dust and Balmer absorp- 
tion worsens this discrepancy. This real vertical offset and 
steeper slope suggests some environmental effect on current 
star formation. 

Using the SFR data from lSkiUman et al.1 (120031 ) for the 
Sculptor Group dwarf galaxies, and a M,/ Lb value of 0.6 
(the average of the _B-band SFR and BdJ values from Ta- 
ble [2] for dl's) to estimate their stellar masses, we find that 
the Sculptor Group dwarfs lie along the same M«-SFR re- 
lation as LG d warfs . In addition, using the Ha data from 
iGavazzi et al.l (119981) for the Goma cluster dwarf galaxies, 
the iKennicutti (| 19981 ) relation between Ha and SFR , and 
the same M^/Lb value to estimate their Af», we find that 
the Coma cluster dwarf galaxies roughly lie on the SFR-M* 
relation of the SDSS (see Fig. I17|l . Since we expect that the 
average SDSS galaxy lies in denser environments than the 
Local and Sculptor Groups, these observations suggest that 
the SFR-Mt relation for dwarf galaxies depends on environ- 
ment. 
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Figure 17. SFR data for the Sculptor Group dwarf galaxies 
(solid red triangles) and the Coma Cluster dwarf galaxies (cyan 
stars) plotted together with the LG dl and Tr gal axies (blue and 
green). Data for the Sculptor Group is take n fromlSkillman et alj 
1120031 ). and for the Coma Cluster, we used lGavazzi et al.l 1119981) . 
We assum e a M^ /Lb value of 0.6, and adopt the prescription of 
iKennicutti 111998 1) to estimate SFR from Ha measurements. The 
solid line is roughly (estimated by eye) the mode line of the dis- 
tribution of SFR in bins of log M« for the SDSS, while the dotted 
lines roughly mark the lines where the condi tional likelihood of 
SFR a given M* is equal to 0.02 (see Fig. 17 of lBrinchmann et al.l 
|2004) . 



6.4 The "Fundamental Line" of the LG Dwarf 
Galaxies 

We have used PCA to quantify the distribution of the LG 
dwarf galaxies in the parameter (log) space of the structural 
quantities alone, and also of the structural + SF quantities. 
We make use of the output of PCA to display the full di- 
mensionality and extent of the distribution of the galaxies 
from the most useful orthogonal angles. 

Using PCA, we find that the LG dwarf galaxies (all 
types) form one linear distribution in structural space 
Mt,/it (or Rt), V. Despite the fact that dE's tend to have 
higher Z than dl's, the distribution of the entire dwarf pop- 
ulation remains linear when we add Z to the parameter 
space. Since this fundamental line resides in the same space 
with consistent slopes as those predicted by supernova feed- 
back models ( M,, /j<. (or R,\ V, Z space), we call this space 
"SNF" space (|Dekel fc Wooll2003l ). 

Since the star formation properties of dE's and dl's dif- 
fer appreciably, especially Z, we performed PCA on dE's and 
dl's separately. We find that the dE's form an even tighter 
linear distribution in SNF space. The dl's are also linear in 
SNF space, and are even tighter in the space spanned by 
Mt,^it,V,Z,M*. 

The linearity of these distributions suggest that one pri- 
mary parameter governs their physics. Furthermore because 
dl's tend to be more massive than the dE's, halo mass may 



be one of the drivers of the distinction between the two 
types. 

Comparison with Previous Work 

PB02 found a "fundamental line" in the parameter 
space (to which we will refer as "PB") of total mass-to- 
light ratio Mtot/Lv, me tallicity [Fe/H], and central surface 
brightness Ev (see also ISimon et al.ll200(j ). For Mtot/Lv, 
PB02 used total masses, estimated dynamically, and total lu- 
minosity. PB space can be derived from SNF space through 
the following relations: 

Mtot/Lv = V^R/L oc V'^R/M, oc V'^/^/alJZ 
Ev oc — 2.51og/x, 



[Fo/H] (X log Z. 



(4) 



PB space is simply SNF space collapsed into three di- 
mensions. We have improved on their analysis in two ways. 
Firstly, we included more dPs in our sample which con- 
tribute enough scatter so that it is no longer linear in PB 
space. Secondly, we have broken PB space into three inde- 
pendent structural parameters (M, V, and R), and metal- 
licity Z. PB space relates the three structural parameters by 
only one relation, namely log Mtot/ii' oc 0.67 log /i* (PB02), 
which translates to one projection of the distribution. A lin- 
ear distribution in the three-space of the structural parame- 
ters should be described by two orthogonal projections which 
we graphically and q uanti tatively present. 

iBernardi et al.l ((2003) studied the fundamental plane 
of SDSS early-type galaxies. They find that these galax- 
ies obey R oc fj~^-^j~o.75^ with slight variations between 
filters, where I = {L/2R~'^) and a is the velocity disper- 
sion. Although the plane equation describing the LG dE's 
R* oc y2i±0-7^o.6±o.2 jg consistent with that of the SDSS 
early-type galaxies, the LG dE's are more accurately a linear 
distribution in the space spanned by 7?*,a, /i* (a permuta- 
tion of M* , V, n, structural space ) . This is consistent with 
the findings of' Bender et al.l (11992') who show the line of dE's 
protruding from the plane of elliptical galaxies in "k" -space, 
whi ch is a permutation of structural s pace. 

IZaritskv et"aLl (|2006al ) (and later IZaritskv et"aLll2007l ) 
find that the LG dE's lie along a "fundamental manifold" 
of spheroidal galaxies in the space spanned by i?e, cr, In and 
M/L. (They note that logAf/L oc log^ a and so add the 
non-power-law term log^ a to their plane equation instead 
of logM/L.) One can easily reduce these four parameters to 
three independent structural parameters (L, a and R or [i) 
through Eq. Q spanning a 3D space in which the LG dE's 
form an even tighter line. 

So in summary, while the dl's seem to extend the struc- 
tural scaling relations of giant late-type galaxies to lower 
masses, the dE's seem to depart from the fundamental plane 
of giant early-type galaxies to form a distribution described 
by a one-parameter family of equations. 



7 CONCLUSION 

Our analysis of structural and spectroscopic data for LG 
galaxies has provided important constraints for modelling 
the physical processes that govern the formation of dwarf 
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Table 9. Scaling Relations 



vs. log A/« 


Sbis 


6bis 


O-bis 


"ID 


bin 


0"1D 


All LG dwarf galaxies 












log M* 


0.61 ± 0.05 


2.94 ± 0.35 


0.38 


0.49 ± 0.05 


3.80 ± 0.35 


0.41 


logK. 


0.33 ± 0.02 


-2.81 ± 0.15 


0.21 


0.28 ± 0.03 


-2.44 ± 0.20 


0.21 


log V 


0.27 ± 0.01 


-0.51 ± 0.11 


0.10 


0.26 ± 0.02 


-0.39 ± 0.12 


0.11 


logZ 


0.40 ± 0.04 


-6.07 ± 0.30 


0.27 


0.29 ± 0.04 


-5.33 ± 0.26 


0.28 


dE galaxies 














log (J* 


0.74 ± 0.05 


2.25 ± 0.37 


0.27 


0.66 ± 0.06 


2.76 ± 0.40 


0.35 


logfl. 


0.26 ± 0.05 


-2.36 ± 0.30 


0.15 


0.17 ± 0.03 


-1.80 ± 0.20 


0.20 


logV 


0.29 ± 0.02 


-0.59 ± 0.17 


0.09 


0.27 ± 0.03 


-0.47 ± 0.20 


0.15 


logZ 


0.40 ± 0.05 


-5.87 ± 0.29 


0.10 


0.37 ± 0.05 


-5.68 ± 0.33 


0.20 


dl galaxies 














logM* 


0.70 ± 0.08 


2.01 ± 0.69 


0.37 


0.51 ± 0.09 


3.52 ± 0.81 


0.53 


logB. 


0.41 ± 0.04 


-3.44 ± 0.41 


0.20 


0.32 ± 0.04 


-2.69 ± 0.39 


0.27 


logV 


0.26 ± 0.02 


-0.43 ± 0.16 


0.08 


0.24 ± 0.02 


-0.29 ± 0.16 


0.14 


logZ 


0.49 ± 0.06 


-7.09 ± 0.49 


0.24 


0.43 ± 0.05 


-6.64 ± 0.42 


0.27 


log M, 


1.18 ± 0.08 


-11.69 ± 0.61 


0.25 


1.14 ± 0.08 


-11.35 ± 0.66 


0.35 



galaxies. In particular, successful models will need to re- 
produce the fundamental line in structural space +Z (i.e. 
SNF space) obeyed by all dwarf galaxies. One such model is 
the sim ple supernova feedback prescription of iDekel fc Wool 
(12003) which quite accurately predicts the slopes of the scal- 
ing relations using M* as the one primary parameter deter- 
mining the distribution of the data. 

What supernova feedback does not deal with is the 
zero-points of the scaling relations, particularly the differ- 
ent zero-points between early and late-type dwarf galaxies. 
Future successful models must include secondary physical 
processes that distinguish between early and late-type dwarf 
galaxies, producing the separate higher-dimensional funda- 
mental lines that they are observed to follow in SNF and 
Ma,, /J.*, V, Z, M-t spaces. The graphical display of PC A out- 
put as we have done in this paper can also help in visualising 
and understanding scatter in model dwarf galaxies. 

Several unexplained challenges include: 

• the different normalisations of the metallicity and sur- 
face density scaling relations for different Hubble type; 

• the presence of recent star formation episodes in 
late type galaxies, and hence the secondary fundamen- 
tal line that dFs follow (i.e. in the parameter space of 
M*,/i*,l/, Z, M.); and 

• the dependence of the dwarf galaxy SFR scaling relation 
with environment, particularly the steepening of the relation 
in lower density environments, and the order of magnitude 
drop in SFR in lower density environments. 
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APPENDIX A: M,/L FROM SSP MODELS 

iMated (|l998l . hereafter M98) has attempted to reconstruct 
the relative star formation rate (SFR) as a function of age 
for many of the dwarf galaxies. We define a function vE'(t) 
which is simply the SFR as a function of t, where t = is the 
time of the first burst of star formation, and call ^(t) the 
"star formation history" or SFH. We normalize '^{t) such 
that its integral over all time is 1. 

iBruzual & GharlotI l|2003l . hereafter BG03) have mod- 
elled the spectral evolution of simple stellar populations 
(SSP), i.e. instantaneous bursts of star formation, using 
evolution a ry tr acks from their "Padova 1994 library" and 
iGhabrierl (120031 ) IMF (which produces very similar results 
for Mf/L as the Kroupa IMF - see their paper). BG03 give 
stellar mass Mt,{t) and absolute V magnitudes as a func- 
tion of time of a stellar population with a mass normalised 
to IMq at the time of the burst, for a population of a 
given metallicity Z. Model results cover the range Z=0.0001, 
0.0004, 0.004, 0.008, 0.02, and 0.05. 



However, 



also 



function of time. 

iBinnev & Merrifieldl l| 19981 ') describe different chemical 
evolution models where Z is roughly proportional to 
Mt/Mgas- To simplify the calculation, we adopt Z ex Af, , 
with 



Z{t) = Zu 



*(f')di' 



(A-1) 



where Ztot is from the measured metallicity (Eq. ((T| in the 
main text). Since Mgas declines with time as the gas be- 
comes stars, Z{t) increases somewhat faster with the time 
of the burst than our simplified Z{t) in Eq. ([A-ip . However, 
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our resulting M,/L ratios are not very diflerent from the 
model with constant Z (the largest difference in log Mt/L 
being 0.05, or less than 1% of log Af«) and we conclude that 
our simplified function Z{t) will not significantly overesti- 
mate M«. 

Given a SFH of a galaxy, ^(f), we can use the BC03 
SSP models to predict the final total luminosity L{to) of the 
galaxy from a convolution of the SFH with the luminosity 
due to V(t,Z(t)): 



Lx = 



^{t')hit',Z{t'))dt' 



(A-2) 



where t = to is the present, and l\{t, Z(t)) is the luminosity 
of the burst for a particular band A in physical solar units 
such that: 



Zv(t,^(i)) = lO°-*[^^'0-^(*'^™l. 
'^{t) is normalised 



(A-3) 



These are simply the Pearson correlation coefficients for the 
m parameters. 

After this reduction of the data, the covariance or corre- 
lation matrices are passed to the PGA. Let the matrix C be 
either Gov or Cor. Then PGA runs on C so as to calculate 
the m X m matrices V and D such that 



CV = VD 



(B-7) 



where V is the product of the Jacobi rotation matrices that 
give the diagonal matrix D. The elements of the diagonal 
matrix Dj,j are thus the eigenvalues of C, and the column 
vectors Vj are the eigenvectors of C. The Jacobi rotations 
are perf ormed by the su broutine jacobi from Numerical 
■Recipes. IPress et al.lll992l . §11.1). Our program for PGA is 
a modified version of a routine kindly provid ed to us by 
Ofer Lahav (described in iMadgwick et al.ll2002l ). Our mod- 
ifications allow for error estimates in the eigenvalues and 
vectors. 



*(t) df = 1 



(A-4) 



so that Lx is the total light in solar units for the star burst. 
Similarly, the total stellar mass is: 

M,,tot= I *(t')A'/*(io-i',^(io-i'))di'. (A-5) 



APPENDIX B: THE MECHANICS OF PCA 

Given a data set with m parameters, the PGA quantifies the 
distribution of the data by performing a series of rotations 
on these original basis vectors of the parameter space. 

Before applying the PGA, the data are reduced in the 
following way: If the original data set contains n galaxies 
and m physical parameters, then we can construct an n x m 
data matrix: 



^i.' 



l..n, j — l..m 



(B-1) 



To simplify our calculations, the matrix is centred about the 
means of each parameter: 



^i.' 



(B-2) 



In order to better estimate the strength of the correlations, 
the matrix can be "standardised" by dividing the elements 
by the standard deviations of the parameters: 



^1, 



(B-3) 



where 



a, = 



ELiC 



ij 



^i) 



^_,^2 



n-1 
Then the covariance or correlation matrix is calculated 



(B-4) 



Gov,,fc = — !-^ Y" yl^A 



8 = 1 



i.j-^i,k 



(B-5) 



and 



B-0.1 Error Estimates in PCA 

The error estimates in PGA were estimated using a "boot- 
strap" method. If the original data set has n galaxies, 
the "bootstrap" method creates A^b new sets of data Afc, 
k = 1,A'^B by randomly selecting n galaxies from the orig- 
inal data set (allowing each galaxy a chance to be selected 
more than once). Then the procedure described above is 
performed on the new data sets A^ to produce Vj, and D^, 
fc = 1, TVs. The bootstrap standard error is 



2 

o"s,v 



1 



A^B 



JVb 



Vfe-Vi3)' 



where V is 



JVb 



Vf 



hT.^" 



Nb 



(B- 



(B-9) 



Analogous relations yield aB,u- 

For our Bootstrap analysis, we used Nb = 1000 to pro- 
vide sufficient statistics for estimating the bootstrap stan- 
dard error. 



B-0.2 Projecting the Data in the New Basis 

The elements of the eigenvectors give the strength of the 
dependence of the vector on the original basis parameters. 
The original data d for the ith galaxy (row) may be trans- 
formed by 



Y = C,V, 



(B-10) 



where the components of Y can be seen as the coordinates in 
the new orthogonal base. Thus for the unstandardised PGA, 
the components of the principal eigenvector give the scaling 
relations between the original parameters. Since these scal- 
ing relations depend on the range of the original parame- 
ters, it is often useful to standardise the data according to 
Eq. (|B-3|I in order to eliminate any bias toward parameters 
with the largest range. 



Cofj.fc = z T y^ '^'j'^t' J' ^ = I--'" 



n- 1 ^-^' 

j=i 



(B-6) 
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